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In supersymmetric quantum mechanics, shape invariance is a sufficient condition for solvability. 
We show that all conventional additive shape invariant superpotentials that are independent of h 
obey two partial differential equations. One of these is equivalent to the one-dimensional Euler 
equation expressing momentum conservation for inviscid fluid flow, and it is closed by the other. 
We solve these equations, generate the set of all conventional shape invariant superpotentials, and 
show that there are no others in this category. We then develop an algorithm for generating all 
additive shape invariant superpotentials including those that depend on h explicitly. 
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Supersymmetric quantum mechanics (SUSYQM) [1- 
3] is a generalization of the ladder operator formalism 
usually attributed to Dirac [4] . It makes use of first order 
differential operators A ± = =p h + W(x, a), where the 
superpotential W(x, a) is a real function of x, and a is a 
parameter. We define two partner hamiltonians 
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H T =A ± A^ = -h 2 — + V T (x,a) , 
ax'' 



(1) 



where partner potentials V±(x,a) are related to the su- 
perpotential by V±(x,a) = W 2 (x,a)±h^^. Part- 
ner hamiltonians [5] have the same energy eigenvalues 
(except for the ground state); i.e., E, : 
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Eq ' = 0. Eigenfunctions of H T are related by: ip. 



n+l 



EV 1 and 



(+) 



A {x,a)ipn ' and A + (x,a)^p^i' 1 oc tp. 
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Thus, if the 

eigenvalues and the eigenfunctions of i2_ are known a 
priori, they are automatically determined for H + as well. 

If the partner potentials V±{x,a) obey the "shape in- 
variance" condition [6, 7], 



V+(x, a ) + g{a ) = V-(x, a%) + g(ai) , 



(2) 



then the spectrum for either Hamiltonian can be derived 
without reference to its partner. This is due to the exis- 
tence of an underlying potential algebra [8-10]. 

Shape invariant partners have the same form except for 
the value of the parameter aj, where a% is a function of do; 
i.e., ai = /(do). The energy eigenvalues of H_(x,ao) are 
given by E^~\ao) = g(a„) - g(a ), where a„ = f n (a ) 
indicates / applied n-times to ao [3]. If the parameters 
differ only by an additive constant: o^+i = a, + h, the 
potentials are called "additive" or "translational" shape 
invariant. All exactly solvable potentials discovered thus 
far that are expressible in terms of known functions are 
additive shape invariant [3, 11]. Several groups found 
these potentials by imposing various ansatzes [10, 12-14]. 

Important correspondences exist between quantum 
mechanics and fluid mechanics [15]. SUSYQM is well 
known to have a deep connection with the KdV equa- 
tion [16-20], a nonlinear equation that describes waves in 



shallow water. We now prove that every additive shape 
invariant superpotential that does not depend on H ex- 
plicitly corresponds to a solution of the Euler equation 
expressing momentum conservation for inviscid fluid flow 
in one spatial dimension. We use this correspondence to 
find a systematic method which (1) yields all known such 
superpotentials for SUSYQM and (2) shows that no oth- 
ers exist. We then extend this method to general additive 
shape invariant superpotentials including those that de- 
pend on h explicitly [14]. 

Writing Eq. (2) in terms of the superpotential yields 



W 2 (x,a ) 



h dW{x,a ) 
dx 

W 2 (x, ai ) - h 



ff(oo) 

dW(x, ai) 
dx 



+ g( ai ) . (3) 



Eq. (3) is a difference-differential equation relating the 
square of the superpotential W and its spatial derivative 
computed at two different parameter values: (x, ao = a) 
and (x,ai = a + h). This equation must hold for any 
value of h. At this point we consider only superpoten- 
tials that do not depend explicitly on h, but only depend 
on H through the parameter a; we will call this class "con- 
ventional" . We will later consider general superpotentials 
that may depend on h explicitly. We show that the shape 
invariance condition (Eq. 2) can be expressed as a local 
non-linear partial differential equation; i.e., all terms can 
be computed at the same point (x,a). This will provide 
a systematic method for finding superpotentials. 

Since Eq. (3) must hold for any value of K, if we ex- 
pand in powers of h, the coefficient of each power must 
separately vanish. Provided that W does not depend 
explicitly on h, this expansion yields 
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da 

W{x,a) = , 



1 dg{a 
2 



da 
n > 3 . 



(4) 
(5) 
(6) 



da n ~ 1 dx 

Thus, all conventional additive shape invariant superpo- 
tentials are solutions of Eqns. (4-6). Although this rep- 
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resents an infinite set, note that if equations at 0(h) and 
0(h 3 ) are satisfied, all others automatically follow. 

Replacing W by —u, x by t, and a by a; in Eq. (4), we 
obtain: 

, d , . du(x,t) Idq(x) 

This is equivalent to the equation for inviscid fluid flow 
in the absence of external forces on the body of the fluid: 

du(x, t) , . „ . . Vp(x, t) 
_Lli +ti (x )t ).V«(x )t ) = -^ (8) 

in one spatial dimension with the correspondence ^ = 

i ^| , where u is the fluid velocity at location x and time 
t, p is the pressure, and p is the local fluid density. Equa- 
tion (8) is one of the fundamental laws of fluid dynamics, 
and was first obtained by Euler in 1755 [21]. Thus, all 
conventional shape invariant superpotentials form a set 
of solutions to the one-dimensional Euler equation. 

I 



Note that Eq. (8) is not closed as written. In fluid 
dynamics this equation is generally supplemented by 
the continuity equation expressing conservation of mass, 
along with an equation of state and/or the energy equa- 
tion and boundary conditions. These additional con- 
straints do not apply in SUSYQM. Instead, Eq. (6) sup- 
plies the additional constraint. 

Equation (6) is satisfied for all n > 3 as long as 

The general solution to Eq. (9) is 

W(x,a) = a-X 1 (x)+X 2 (x)+u(a) . (10) 

Substituting this into Eq. (4), and collecting and label- 
ing terms based on their dependence on X\ and X2, we 
obtain 



^ +(- d ~t )+ a + H^) + ^ + •§)*-*<.>• <») 

Term#l v ' Tcrm#3 v ' ^~^~^ K v ' 

Term#2 fr Term#4 Term#5 Term#6 



where H{a) = -u$± - ± To find all possible solu- 
tions, we begin by considering special cases of Eq. (11) 
where one or more of the terms X^x), X 2 (x), or u is 
zero. After considering these cases, we will show that 
all solutions can be reduced to one of these cases. In 
our nomenclature, lower case Greek letters denote a- and 
cc-independent constants. 

Case 1: X 2 and u are not constants, X\ is 
constant. In this case, let X\ = p. Then W = 
pa + u(a) + X 2 {x). If we define u = u(a) + pa, we get 
W = u + X 2 . So this case is equivalent to Xi = 0. Then 
terms 1, 3, 4, and 6 each becomes zero, and Eq. (11) 
becomes + ^X 2 = H(a). Since X 2 must be in- 

dependent of a, and H(a) must be constants. This 
yields u = aa + (3 and — + aX 2 = 8. The solution is 
X 2 (x) = f+rye^. Therefore, W = aa + f3 + f + 77 e ax . 
Defining a = —1, we obtain W = A — Be~ x , where 
A = f3 — a — 9. This is the Morse superpotential. 

Case 2: X\ and u are not constants, X 2 is con- 
stant. Following a similar procedure, this case is equiv- 
alent to X 2 = 0. Depending on the values of constants 
of integration, this equation yields the Rosen-Morse I, 
Rosen-Morse II, Eckart, and Coulomb superpotentials. 

Case 3: X\ and X 2 are not constants, u = pa + v. 
We define X\ = X\ + p and X 2 = X 2 + v, making this 
case equivalent to u — 0. Depending on the constants of 
integration, this yields the Scarf I, Scarf II, 3-D oscillator, 



and generalized Poschl- Teller superpotentials. 

Case 4: X 2 is not constant, X\ and u are con- 
stant. If X\ ^ we get Morse, and X\ = generates 
the one-dimensional harmonic oscillator. 

Case 5: X\ is not constant, X 2 and u are con- 
stant. This yields special cases of Scarf I and Scarf II, 
and the centrifugal term of the Coulomb and 3-D oscil- 
lator. 

Case 6: X\ is constant, X 2 is constant. In this 
case, the superpotential has no x-dependence, regardless 
of the value of u. This is a trivial solution corresponding 
to a flat potential, and we disregard it. 

These special cases generate all known conventional 
additive shape-invariant superpotentials [3, 11], as shown 
in Table I. 

Now that we have considered these special cases, we 
can systematically obtain all possible solutions. H (a) is 
independent of x. Therefore, when any solution is substi- 
tuted into Eq. (11), it must yield an x- independent sum 
of terms 1-6. There are many ways in which these terms 
could add to a sum independent of x. We begin with the 
simplest possibility, in which each term is individually 
independent of x. In this case, term 3 states that X\ 
must be a constant, independent of x. In addition, term 
1 dictates X\X 2 must be constant as well. These two 
statements can only be true if X 2 and X\ are constant 
separately; this reduces to the trivial solution of case 6. 
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Therefore, assuming that each term is separately in- 
dependent of x yields only the trivial solution. How- 
ever, there is also the possibility that some of the terms 
depend on x, but when added to other terms, the in- 
dependence cancels to yield a sum that is independent of 
x. If a group of n-terms taken together produces an x- 
indcpcndcnt sum, and if no smaller subset of these terms 
add up to a sum independent of x, we call this group 
"irreducibly independent of a;." 



Name 


superpotential 


Special 
Cases 


Harmonic Oscillator 




Xi = u = 


Coulomb 


e A e+i 

2(1 + 1) r 


X 2 = 


3-D oscillator 


V- — 

2 r 


u = 


Morse 


A - Be~ x 


Xi =0 


Rosen-Morse I 


-Acotx - -f 


X 2 = 


Rosen-Morse II 


A tanh x + -f 


X 2 = 


Eckart 


-Acothx + -§ 


X 2 = 


Scarf I 


A tan x — Bsec x 


u = 


Scarf II 


A tanh x + Bsech x 


u = 


Gen. Poschl-Teller 


A coth x — Bcosech x 


u = 



TABLE I. The complete family of conventional additive 
shape-invariant superpotentials. 



If, for example, term 2 depends on x and term 5 
depends on x, but the sum of these two terms is x- 
independent, then we consider the set of terms {2, 5} to 
be a two-term set that is irreducibly independent of x. 
Let us investigate this example further. 

In this case, — + ^X 2 is independent of x. How- 
ever, Xi and must each depend on i, or this would 
be reducible. Since term 2 does not depend on a, jj^ 
must be constant for the x-dependence of terms 2 and 
5 to cancel. So u = 5±a + 5%. Substituting this into 
terms 2 and 5 yields — + SiX 2 = 5 3 . The solution is 
X2 = |^ — Sie SlX . For this solution to work, the sum of 
the remaining terms 1, 3, 4, and 6 must also be indepen- 
dent of x. We first ask if this could be true by making 
all of the remaining terms each independent of x. This 
is only possible if X\ = 0. Thus, the combination where 
terms 1, 3, 4, and 6 are each individually independent 
of x and {2,5} is an irreducibly independent set is an 
example of Case 1 above (since X\ = 0) and yields the 
Morse superpotential. 

We continue in this manner, checking whether each 
two-term irreducible set yields solutions when combined 
with the remaining terms each independent of x as in 
the example above. In each case, we find either that the 
equation reduces to one of the special cases examined 
earlier, or that no solution is allowed (for instance, term 



1 and term 3 cannot be irreducibly independent of x since 
one is independent of a and the other is linear in a). 

Once these combinations are exhausted, we consider 
combinations of two-term irreducible sets with other two- 
term irreducible sets as well as single-term constants. 
Then we examine three-term irreducible sets, all the way 
up to the full six-term equation. 

As a final example, we check whether there are any so- 
lutions for the full six-term irreducible set. We note that 
the first two terms are independent of a, while terms 
3 and 4 are linear in a. We do not know a priori the 
functional form of u. However, we do know that any in- 
dependence in terms 1 and 2 cannot be canceled by terms 
3 and 4, since terms linear in a cannot cancel terms inde- 
pendent of a. For an irreducible set, the sum of the first 
two terms must have an x-dependence that is canceled 
by a-independent terms from u and ^ in terms 5 and 
6, and terms 3 and 4 must have an x-dependence that is 
cancelled by terms linear in a. Since term 5 contains 4 s , 
it could include terms independent of a, terms linear in a, 
and/or other forms of a-dependence. However, it cannot 
fully cancel the rr-dcpcndcnce of the first four terms or 
the set will be reducible. 

We conclude that the only way for the solution to be 
irreducible is if u + = + fia + v for constants 
H, v, and e This gives" u = f + (f + f) + By 
collecting terms of the same power in a, we find that the 
terms proportional to ( a l v y mandate that 7 (X2 + vX\) 
is a constant. This leaves only two possibilities. First, 
7 = 0, in which case u depends only linearly on a, and 
this reduces to Case 3 above. Otherwise, X2 + vX\ must 
be a constant. In this case, X\ differs from X 2 by only 
a multiplicative constant; by shifting the zero of a, we 
can absorb X2 into Xi, and this reduces to Case 2 above. 
Thus, any solution for the full set of terms can be reduced 
to one of the special cases. 

By examining all possible combinations of terms, we 
have found that no new solutions are admitted by any 
combinations that are not included as one of our special 
cases. Thus, we have found all known additive shape- 
invariant superpotentials that do not depend explicitly 
on h, and have proven that no more can exist. 

However, a new family of "extended" shape-invariant 
potentials was recently discovered by Quesne [14], and 
expanded elsewhere [22]. These potentials are generated 
from our system by generalizing our formalism to include 
superpotentials that contain h explicitly. In this case, we 
expand the superpotential W in powers of h: 

00 

W{x, a,h) = J2 h n W n (x, a) . (12) 

n=0 

Substituting Eq. (12) in Eq. (3), significant algebraic 
manipulation yields 
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E ft " 



71=1 



E^- + ^-EE 



Tl-l 



1 fln-s " * 1 



(n- s)\ da r < 



(k — 1)! da k dx 



Qk+l / 1 gn g - 



0. 



Lfc=o s=0 k=o v " fc=0 

As this must hold for any value of h, the following equation must hold separately for each positive integer value of n: 
<)\V n ■ v-x-- 1 0" - ^ I d k 

W k W n - k + 

k=0 



dx 



n s 1 
E E ( n _ 



s=0 fc=0 



(n - s)! da"-" 



W fe VF s _ fc 



n 1 

^ (fc-1)! da k ~ x dx 



W„_ fe 



i , N 

— : 7T — 1=0. 13 



For 



1, we obtain 



dx da 



(W^+g) = 0, 



(14) 



yielding 2& = ^ (W 2 + g) for > 1. We have 
shown that all conventional superpotentials W = Wo are 
solutions of this equation. Higher order terms can be 
generated from applying Eq. (13) for all n > 1. 

As an example, we choose the 3-D oscillator solution: 



W 



For n = 2, the expansion yields 



dWi _ d_ 
dx da 
and for n = 3, we obtain 

dW 2 d (2W W 2 + W?) 



dx 



da 



{WoWx) = 



1 d 2 WoW x 2 d 3 W a 

2 da 2 + 3 da 2 dx ~ ' 



These two coupled equations are solved by W\ = and 
W2 = (4xw)/(2a + x 2 u;) 2 . The next order equations are 
solved by W 3 = and VF 4 = (4xw)/(2a + x 2 cj) 4 . Gener- 
alizing these, we get 



Wo 



1 



W 2n+1 = 0; W 2n = (42xj)/(2a + :r 2 u;) 



2n 



yielding a sum that converges to 



2ujxh 



2uoxh 



1 a 
VF (x, a, ft) = -cjx hi , 

2 a; \ljx 2 + 2a - h ux 2 + 2a + h 

With the identification a = (£ + l)h, and ft = 1, 



W 



LUX 

~2 



£ + 1 



+ 



2ux 



2ujx 



ujx 2 + 2£+l lux 2 + 2£+3 



This is the extended superpotential found by Quesne [14]. 

We have thus obtained a system of partial differential 
equations that must be satisfied for all shape-invariant 
superpotentials. For conventional cases that do not de- 
pend on h, we have shown that the shape invariance con- 
dition is equivalent to an Euler equation expressing mo- 
mentum conservation for fluids and an equation of con- 
straint. For extended cases in which the superpotential 
depends explicitly on h, we developed an algorithm that 
is satisfied by all additive shape invariant superpoten- 
tials. 
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